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1 Introduction 

According to AdS/CFT correspondence, any global symmetry at the boundary theory is 
lifted to a local symmetry in the bulk [1, 2]. The gauge symmetry is essential to reduce 
the degree of freedom which is enlarged by going into one higher dimension. The phys¬ 
ical goal in holography is the boundary quantities which does not know the presence of 
higher dimension or gauge degrees of freedom, while we use the tools in the bulk theory. 
Therefore the gauge invariance of a physical quantity is a critical issue for the validity 
of the AdS/CFT. Also tracing the gauge invariance gives much intuition on the way how 
holography actually works, especially how global symmetry is encoded in the local gauge 
symmetry. 

One can find gauge invariant combinations of the fields, and express the physical 
quantities in terms of such master variables, however, it is not always easy to find such 
gauge invariant combination. Even in the case they are available, it is not very convenient 
to use such fields, especially if many fields are coupled, because the physical quantities are 
defined in terms of the field variables which are formally gauge dependent. For example 
[2], energy momentum tensor and chemical potential is defined in terms of metric/gauge 
field which are not gauge invariant. Similarly, heat currents can be related to the metric 
perturbation defined only in a specific gauge where time period has definite relation with 
temperature. 

In recent works [3, 4], based on [5, 6], we developed a systematic method to numerically 
calculate the Green’s functions and all AC transports quantities simultaneously for the 
case where many fields are coupled and there are constraints due to gauge symmetry. 
Although we have tested the validity of the procedure by showing the agreement of zero 
frequency limits of AC conductivities with the known analytic DC conductivities [7-9] we 
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still think that we need to prove the gauge invariance of our procedure as a matter of 
principle. We found that the bulk gauge invariance is intimately related to the holographic 
renormalisation. Although the local counter terms were introduced to kill the divergences, 
it also kills most of gauge dependence. 

Furthermore, there is a residual gauge symmetry(RGS) even after we fix the axial gauge 
9rx = 0. While equations of motion can be written in terms of the the gauge invariant 
master fields fP/,, CP X (3.8), it turns out that the quadratic on-shell action, the generating 
function for two point retarded Green’s functions, can not be written as such. However, 
we prove that the Green’s functions are still invariant under such a symmetry. 

There is a mismatch in the degrees of freedom in the bulk and those at the boundary: 
there are only two independent bulk solutions satisfying the in-falling boundary conditions 
while we need three solutions at the boundary since there are three independent source 
fields. The RGS is the one that resolves the problem: since it cannot satisfy a proper 
boundary condition, it is not a proper gauge symmetry but a ‘solution generating sym¬ 
metry’. It generate the desired solution at the boundary and therefore we should accept 
its bulk counter part as a new physical degree of freedom as well although it can not sat¬ 
isfy the infalling boundary condition (BC). By extending the RGS such that it satisfies 
infalling boundary condition at the horizon, we can make the bulk solution more natural 
in the sense that it satisfies the infalling BC. With such solution we can also understand 
the problem in the context of general structure of holography, namely the correspondence 
between a global symmetry at the boundary and the local gauge symmetry in the bulk. 


2 Action and background solution 

Let us first briefly review the system we will discuss, which has been analysed in detail in 
[3, 7, 10]. The holographically renormalised action(S ren ) is given by 

Sren = Sem + + S c , (2.1) 


where 


Sem = / d A Xy/^g 
J M 


R-2A- -F 2 

4 


- 2 / d 3 xV^/K , 

JdM 


( 2 . 2 ) 


is the usual action for charged black hole in AdS space(A < 0) with the Gibbons-Hawking 
term and 


S 4’= [ dix V~9 

■I M [ Z 


(2.3) 


is the action for two free massless scalars added for a momentum relaxation effect. S c is 
the counter term 


S c = Vc 



dx 3 \G-7 


(- 4-RM + 


\ X 


(2.4) 
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which is included to cancel the divergence in Sem + <SV Here we introduced r) c to keep 
track of the effect of the counter term. At the end of the computation we will set r] c = 1. 
The action (2.1) yields general equations of motion 1 


Rmn = t;9mn 


R- 2A - -F 2 - - 
4 2 


2 \ i i 

E(W 2 ) + - 2 E dM^idNipi + 2 F m P Fnp , (2.5) 


i=i 


v m f mn = 0 , v 2 ^/ = o, 

which admit the following solutions 


dr 2 

ds 2 = G MN dx M dx N = — f(r)dt 2 + —— + r 2 5ijdx l dx^ , 

f{r) 

d2 


/3 2 m 0 , V 2 r l 






P 2 


4r 2 2r 2 


A = fi ^1 — —^ dt 

i'i = Pux 1 = Pdnx l 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


These are reduced to AdS-Reissner-Nordstrom(AdS-RN) black brane solutions when /3 = 0. 
Here we have taken special (3u, which satisfies ^ Y^i=i Pi' Pi = P 2 f° r general cases. 

The solutions (2.7) - (2.10) are characterised by three parameters: ro, 9 , and /3. ro is 
the black brane horizon position(/(ro) = 0) and can be replaced by temperature T for the 
dual field theory: 


/'(ro) 

47T 



/r 2 + 2/3 2 \ 
4 r 0 ) 


( 2 . 11 ) 


Non-vanishing components of energy-momentum tensor and charge density read 


(T tt )=2m 0 , (T xx ) = (T yy ) = m 0 , (J 4 ) = //r 0 . (2.12) 


(' T tt ) = 2 (T xx ) implies that charge carriers are still of massless character. From here we 
set ro = 1 not to clutter. 

1 Index convention: M, N, ■ ■ ■ =0,1,2, r, and - ■ ■ =0,1,2, and i,j , • • ■ = 1, 2. 
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3 Gauge fixing and residual gauge transformation 


To study electric, thermoelectric, and thermal conductivities we introduce small fluctua¬ 


tions around the background (2.7) - (2.10) 

/ °o J. . 

— e-™a x (u>,r), (3.1) 

-OO 

/ oo T. 

—e~ ru}t r 2 ht x (uj,r), (3.2) 

-OO 

/ OO H 

— e~ iwt r 2 h rx (u , r), (3.3) 

-OO 

/ OO J. . 

— e-^xi^r). (3.4) 


The fluctuations are chosen to be independent of x and y. This is allowed since all the 
background fields appearing in the equations of motion turn out to be independent of x 
and y. The gauge field fluctuation^^(f, r)) sources meAiic(Sgt x (t. r), 5g rx (t, r)) and scalar 
field(<5^i(t, r)) fluctuation and vice versa and all the other fluctuations are decoupled. We 
will work in momentum space and ht x (ui,r) and h rx (u,r) is defined so that it goes to 
constant as r goes to infinity. 

By linearising the full equation of motion, we get four equations. However one of them 
can be obtained by the others. Thus we may consider following three equations: 


(V - Ph„) - irih - ir V'4 + ^ ) = o. 


P r 2 f(r) 


Pf(r) 


a'x( r )f( r ) , w 2 a x (r) n(ti tx + iuh rx ) 


+ 


+ 


= 0, 
\ 2/\7 


a "(r) H- 

^ J f(r) f{r) 2 ' /(r 

f(r)f(r)(x'{r) - (3h rx ) + f(r) 2 (x' ~ Ph rx )' + 2 ^ ^ 

r 

+ u 2 x(r) - if3uh tx (r) 


= 0. 


(3.5) 

(3.6) 

(3.7) 


If we differentiate the third equation with respect to r, all equations can be written in 
terms of three variables, CP X ,CPand a x , where 

y x = X r ~ PK X , CP h = h' tx + iuh rx . (3.8) 

Therefore, h rx is a non-dynanrical degree of freedom. Indeed, CP X ,CP h, and a x are invariant 
under a diffeomorphism generated by = (0, C(?’)e _jaj4 ,0,0), under which the fields are 
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transformed as follows: 


Skrx = + V.e r ) = C(r)e ~ iujt , 

(3.9) 

Shtx = ^(V t 6 + V x &) = -iu;C(r)e- iujt , 

(3.10) 

5 X = K{r)e~ iut , 

(3.11) 

5a x = 0. 

(3.12) 


Using this gauge degree of freedom, one may set h rx = 0, which is so called the axial 
gauge. The numerical calculation in [3] has been performed in this gauge. A question arises 
whether the resulting physical quantities are independent of such gauge fixing condition. 

Furthermore, even after we fix h rx = 0, one can still find a residual gauge transfor¬ 
mation which is given by constant £ [11], This residual diffeomorphism doesn’t change 
the gauge fixing condition h rx = 0 and generates constant shift on ht x and y, because the 
equations of motion contain only derivatives of ht x and x and the linear combination of 
them, w\'(r) — i/3ht x (r), which is invariant under 


htx -t h tx + h 0 , and x 


, .0, 
X+i~h 0l 
UJ 


(3.13) 


where ho is a constant. Thus there is one parameter constant solution given by 


cl x — 0 ? 


htx — ho-, x — ^ ho ■ 

UJ 


(3.14) 


which does not satisfy in-falling boundary condition so it is not a physical degree of free¬ 
dom 2 . We call it the residual gauge symmetry(RGS) because it is generated by the zero 
mode of a diffeomorphism generator. This kind of solution was first introduced in [12] 
Why should there be such a residual degree of freedom? It can be traced to the 
difference of the differential equation near horizon and those near boundary. Near the 
black hole horizon (r —> 1) the solutions are expanded as 


hx = (r- l) v±+1 (4? + 4? } (r - 1) + ■ ■ ■) , 

a x = {r- l) u± (a^ + 4 /J) (r - 1) + ■ ■ ■) , (3.15) 

X = (r - l) u± (x (/) + X (II \r - 1) + • • ■) , 

where v± = ±z4w/(—12 + 2/3 2 + /r 2 ) = =piw/(47rT) and the incoming boundary condition 
corresponds to v = By inserting these to the equations of motion, one can easily find 
a linear relations between the zero-th modes: 


[y + 1 ) 4 ? + + Px {1) = 0 , 


(3.16) 


Notice that all other modes are generated by these. Thus there is a well defined constraint 

2 It is a regular solution at future horizon. 
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equation which reduces the degrees of freedom. 

On the other hand, by inserting the expansion near the boundary (r —> oo) 

h - h {0) -l —h (2) 4 . Jl/,( 3) i 

'Hx — Il-tx "T 2 rl tx "T" r 3 'ha- ~r ' ' ' J 

= oi°' ) + —‘ ‘ ‘ , (3-17) 

x=x (o >+ 4 x ( 2 ) +V 3) +" •. 

fy* Aj fy* O 

to the equations of motion, we can not get any relation between the zero-th modes a^\ , 

and ; all °f which are related to the higher modes. More explicitly, 

w(u;x (0) - i?b$) - 2 X (2) = 0, i/3( WX (0) - ^4?) - 2 h$ = 0, (3.18) 

which are evolution equations in r-direction. Therefore, there is no constraint equation. 
Then there is a crisis of mismatch of degrees of freedom and this crisis is resolved by the 
effective residual degree of freedom described above. However, this residual gauge degree of 
freedom raises another issue of invariance of physics under this symmetry. We will address 
this issue at the end of section 5. 

4 Holographic renormalization and gauge invariance 

Now we come back to the question whether physical quantities are independent of the choice 

of the gauge condition h rx (r ) = 0. We will show this by proving that the generating function 

of physical quantities, the on-shell action, is invariant even in the case with h rx (r ) ^ 0. 

(2) 

The on-shell renormalised action to quadratic order in fluctuation fields, S re n , is 


= lim / d 3 x 


Hi (^PfSgrx - )yfr 2 H'i \ + - ^ftiA x 5A' x 


- Htx 




J 

2 r 2 ' 


(4.1) 


+ vc Hi 


r 2 Hi _ PHtx 
2\/7 2V7 


PSipidgtx ( 2 

+ -^7T - \7f 1 Sg ‘■ 


where f(r) = r 2 — ^-We dropped the boundary contribution from the 

horizon as a prescription for the retarded Green function [13]. 3 Near boundary r —> oo, 

3 In fact, the contribution of the incoming solution at the horizon is zero in (4.1), which is real. However, 
for a generating function of retarded Green’s functions, we will take only part of (4.1) as explained below 
(4.3), which is complex. In this case, it turns out that the contribution from the horizon is pure imaginary. 
From this perspective, we should drop the contribution from the horizon. 
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the fluctuation fields in momentum space, (3.1) - (3.4), may be expanded as 


00 


00 h^(/A 
u t \ _ nrx 1 W 1 
h r x{u,r) - 2_^ 


n =0 


n =0 


a 


00 (n)/ \ 

( \ ax ( w ) 
:( w > r ) = 2 ^^r~ 


X(w,r) = ^ 


X (n) M 


(4.2) 


n=0 n=0 

and using the equations of motion, we can obtain a quadratic action as follows 



V 2 / , °°dw 
2 Jo 2vr 


- M 0 ) 4 ? - m4?4 0) 


2m 0 4° ) 4° ) + 


+ (*<°> + (3x< 3 >+/3/,<S) 

+ (ic - 1)| - a 3 (aaS'aSS 1 ) - A 2 )4A| 3, A™ + 4iS£’ji<J>w) 

+ A (*tfx (0) w - 2 ih^h^u + /3 2 4?4!?) 

+ A (-41^4* w “ 4 4® ■4? + *0X (O) 4!? W ~ X (0) X (0) w 2 ) 

- 2m 0 4?4° ) - - 2iw4x4x + /3 2 4M } + iPuhQx^ 


~ - 4/>M + i/3wx (1) 4? 


- ^ 2 x (1) x (0) 




" 

} 

+ 

c.c 


(4.3) 


where the argument of the fields 4 is ui. V 2 denotes volume in x-y space and [c.c] means 
the complex conjugated terms. From here, we will drop the [c.c] term since we want to 
compute retarded Green’s functions [13]. 

The second line is proportional to a gauge invariant combination under (3.13). Fur¬ 
thermore, one of the equation of motion including fi'r) is 


hix ~ pJZ7^2 ( 3iujh tx ~ i/^4 0) ~ 3 /3x (3) ) = 0. 

One can show that (4.4) is equivalent to a Ward identity 


V m (TH + F x v {J x ) - (0 7 )c)>/ = 0, 


(4.4) 


(4.5) 


by using the boundary metric and the other fields in the linear approximation given as 

4 a*°^(w) = ai°\—uj) = ai°\cu) by the reality condition of 5A X . The same notation and reality condition 

apply to all the other fields. 
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follows: 


ds 2 = rt^dx^dx" + 2h^e~ iut dtdx , (T^) = 

F = -iua^e-^dt A dx , ( J^) = <( J (0) ^ + (J (1)/i ) = (/q 0, 0) + (o, a ,W - , o) e~ iu}t 

i’l = (px,py) , (O 1 ) = (o (1)/ ) = (3 X (3) + /3hg>,o) e~ iujt , (4.6) 

where 



m 0 





= ( -2 m 0 hfj - 3 hfj + iu 



—iu)t 


(4.7) 


One may ask why Ward identity of the boundary theory is included in the bulk equation 
of motion. It is not accidental: The translation, x —> x + £o at the boundary theory is 
imbedded into the bulk diffeomorphism x —> x+£(x), which induces the field transformation 
$ —> 4> + <5g4>, which in turn is a special case of general variation, 4? —> 4> + 54). Now the 
equation of motion is coming from the invariance of bulk action SSb = 0 under the general 
variation, while the Ward identity is the requirement of the boundary action under the 
translation 6 ^ 0 Sb = 0. Because AdS/CFT request Sr = Si, at the onshell, the latter is 
contained in the huge tower of equation of motion as a tiny piece. 

The terms proportional to (r/ c — 1) in (4.3) include the divergent terms with A, a 
regularisation parameter, and finite terms without A. A remarkable fact is that with the 
counter term of weight r/ c = 1, not only the divergent terms are cancelled, but also all 
the h rx dependent finite terms disappears from the on-shell action, as we claimed in the 
beginning of this section. 


5 Gauge invariance under the residual gauge transformation 

Our starting point is the action 5 

- - 2- 3+ 3y (0) x (3) 


(2) V2 f 00 du 


2 ./o 2vr 


+ 


(Vx (0) + «A ( ° } ) 


+ c.c, 


(5.1) 


which is still dependent on residual gauge (3.13) even after we set h rx = 0. Since it is 
just a constant shift of the solution 4>, its effects are only shifts of zero-th modes and 
4> , (r) and all of its modes, especially (a^ x \h[^},x^) ■= 11“ are intact. Notice that the 
recurrence relations derived from equations of motion relate higher modes with the zero-th 
modes J a = (a^\ \ x^)- However, all dependences of higher modes on zeroth modes 

is through the gauge invariant combination wx^ — . See, for example, (3.18). Thus 

5 It comes from (4.1) before we get eq. (4.3), for which we have to use the equations of motion. 
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all higher modes are gauge invariant, which makes the gauge invariance of the &(r) intact 
in spite of the complicated dependence of higher modes on the zeroth modes. 

The residual gauge dependence of (5.1) can be understood as follows. The full on shell 
action should be invariant under the residual gauge transformation. However, what we are 
looking at is the quadratic part of the action n, which generates the 2-point function, in 
the expansion of 

SVenl^d*] = SW + sglm + S® [54>] + • • • , (5.2) 

where = ( 6 $^, <54^, <54>/) collectively denotes the sources of the dual field theory, which 
are boundary values of SA^ and 5ipi. Sren^d)] and Srcn^d*] are given as follows: 

sillm = J d 3 x (t (0) ^) + 8^ (jW") + <54> 7 (o^)) , (5.3) 

SW[ 8 $\ = J d 3 x + 8 ^ (+ 8 $! (o«')) • (5.4) 

Under the residual gauge transformation 6 with ho = —iuC ,Oj the variations of these actions 
are 

^Ven[<$$] = V 2 J ^ {Co + 2iu}m 0 h^ + c.cj , (5.5) 

8 S®m = -SS^m + U 2 f ^ {Co (3 /?x (3) - 3 iuh® + + (/ 3 2 - cu 2 ) /$) + c.c} 

(5.6) 

Thus the total variation is proportional to the Ward identity (4.4). Notice that 5 ren is gauge 
( 2 ) 

invariant but 5 r ' e n, which is starting point to derive the Green function, is not invariant by 
itself. Nevertheless physical observables derived from Sr e n are invariant because the Green 
functions are second derivatives of the full on shell action at the zero source limit. 

At this point one can discuss a puzzle in counting degrees of freedom. There are only 
two independent bulk solutions satisfying the in-falling boundary conditions,' while we need 
three solutions at the boundary since there are three independent source fields. Therefore, 
there is a crisis of mismatch of degrees of freedom between the bulk and boundary. What 
solves the problem is the RGS (3.14). We call it RGS because it is generated by the zero 
mode of a diffeomorphism generator. On the other hand, to be a proper gauge degree 
of freedom in the bulk, the diffeomorphism generator should satisfy the proper boundary 
conditions: infalling at horizon and Dirichlett at boundary. The residual gauge symmetry 
generator is a global shift and therefore it can satisfy neither of them. So such a shift by the 
diffeomorphism zero mode is not a true gauge symmetry, while it is a symmetry of the bulk 
equations of motion. In other words, the RGS is a “solution generating symmetry” rather 

'’This transformation changes the sources of the action, 5$^, 5&i . One should note that there are 

non-vanishing transformations for 54>oo and <5<f>o. 

'We have two second order differential equations and one first order one in three variables: a x ,ht x ,X- 
Therefore, there are 5 boundary conditions to fix. If we fix the in-falling boundary conditions for all three 
variables, we are left with two degrees of freedom. We recall the equations (3.15) and (3.16). 
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than a gauge symmetry. Therefore, the gauge orbit of RGS can provide us the necessary 
degree of freedom (d.o.f) near boundary. To match the d.o.f, we need to accept its bulk 
orbit as physical configuration inspite of the fact that the resulting bulk solution does not 
satisfy the infalling BC. 8 One can give a more natural bulk solution by extending RGS to 
a diffeomorphism which satisfies the infalling boundary condition and it is reduced to our 
previous RGS near the boundary. It is generated by = (0, C t {r)e ~‘ lut , 0, 0), with 9 

C {r) = e{f{r)/r 2 )-^ T \ (5.7) 

where / is the metric factor given in Eq. (2.8) and e is a constant parameter. Notice 
that the RGS is the case where £( r ) is constant. We will call this “boundary shifting 
diffeomorphism” (BSD). Now we can understand the degree of freedom mismatch as follows: 
Since it is not satisfying the Dirichlett be, it is still not a proper gauge transformation. 
Notice also that under (5.7), the gauge slice is shifted and some of the gauge fields become 
singular. For the discussion on the treating these issues, we refer the reader to page 24 of 
ref. [9]. 10 This is the reason why the BSD can generate a new solution in the boundary. 
It is precisely the same logic why RGS generate new solution. * 11 Since RGS and BSD 
shift the boundary values of fields, they generate the Ward identity for the translation 
invariance. This is a typical example how a global symmetry is encoded in a local gauge 
transformation and how the apparent paradox of the degree of freedom can be resolved 
because of the holographic correspondence. 12 

6 Basis independence 

In [3], we constructed a formalism to perform the AC conductivities for the case where 
multiple fields are coupled together. We had to choose a basis of initial conditions and one 
can ask whether different choices of basis give the same result. Answering this question 
will also provide an alternative reasoning of gauge invariance. To provide the setup, let us 

8 So far we discussed the degree of freedom mismatch using the RGS, since our formalism in [3] to 
calculate the conductivity is based on it. 

9 We thank the anonymous referee for suggesting to consider this. 

10 It is very tempting to consider BSD as a gauge transformation at least from bulk point of view. If we 
do it, we get to the problem: Its orbit in the boundary generate physical configuration while it does not in 
the bulk, so that crisis of d.o.f becomes real! 

11 This argument is further justified if we consider the numerical calculation starting from the boundary 
instead from horizon. After choosing 3 fields’s values, we can adjust two “expectation values” such that 
we can get infalling boundary values at the horizon. It is easy to show that only when we start from a 
subspace of codimension 1, we get three infalling solution near the horizon. If we start from a point off this 
plane, we get one infalling and two fields which are mixture of infalling and a constant. In this calculation 
the gauge condition h rx = 0 is intact. This demonstrates that we can not impose infalling be for all fields 
at hands. If we do the same numerical experiment for BSD, the picture is following. The BSD generate 
the the orbit and it also move the gauge slice. Now in this case even in the case we start from the off the 
plane, we can get three infalling fields at the horizon. We need to calculate the r-evolution at each ‘gauge 
fixing’ plane which pass through the initial data. 

12 The apparent ‘mismatch’ is due to the difference in viewing the gauge orbit of BSD (or RGS) between 
the bulk and boundary. In the bulk, one could view it as gauge orbit. On the other hand, from the boundary 
theory point of view, there is no gauge structure and the orbit of translation symmetry is physical degree 
of freedom. 
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consider N fields <F a (x, r), (a = 1, 2, • • • ,1V), 

/ ,]dl, 

( 6 . 1 ) 

where the index a may include components of higher spin fields. For convenience, r p 
is multiplied such that the solution d > a (k,r) goes to constant at boundary. In our case, 
(d* 1 , d> 2 , d> 3 ) = (a x , htx, x) and P = 0 for d* 1 , d> 3 and p = 2 for d> 2 . 

Near horizon(r = 1), solutions can be expanded as 

$“(£:, r) = (r-iy a± + <Pi(r -!) + ■■■)> (6.2) 


where a new subscript i is introduced to denote the solutions corresponding to a specific 
independent set of initial conditions. For example, may be chosen as 


Pi — (i> — (A + £0/(i + u )i i) > 
P2 = (l, ~(p - ^)/(l + v), -l) T , 


(6.3) 


where we used (3.16) and v = — zcj/(47tT) as shown below (3.15) for incoming boundary 
condition to compute the retarded Green’s function [13]. Due to incoming boundary con¬ 
dition, (p? determines through horizon-regularity condition so that we can determine 
the solution completely. Each initial value vector <p l yields a solution, denoted by dh(r), 
which is expanded as 

O a 

(A;, r) —> §“ + • • ■ H—/ + • • • (near boundary), (6.4) 

j'Oa 

where S“ are the sozzrces(leading terms) of z-tli solution and 0“ are the operator expectation 
values corresponding to sources(5 a > 1). 

Notice that we have only two solutions while we have a three dimensional vector space 
3 of boundary values J“, a = 1,2,3. To fix such mismatch of degree of freedom, we 
introduce a constant solution 4>o( r ) = §o = (0, l,i(3/u) J along the gauge-orbit direction of 
the residual gauge transformation so that §“,§2,§o f° rm a basis of 3 • Now § and Q are 
generic regular matrices of order 3. 

The general solution is a linear combination of them: let 


$“(fc,r) = $?(k,ry, 


(6.5) 


with real constants c*’s. We can choose d such that the combined source term matches the 
boundary value J a : 

J a = §“c 4 , (6.6) 

which yields 

ID 

<h a (fc, r) = 4 > “(/c, r)d — > J a + ■ ■ ■ -|—- + • ■ ■ , (near boundary) (6.7) 

j'Oa 
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where, with (6.4) and (6.6) 


n a = ©“c* = ©“(S ~ l ){J b =: C%J h . 


( 6 . 8 ) 


Notice that both 11“ and C b are invariant under the transformation J b —> J b + eSg because 
C b S § = ©“(S- 1 )^ = ©o = 0, where ©g = 0 since it is the sub-leading term of the 
constant solutions. 

A general on-shell quadratic action in momentum space has the form of 



d d k 
(2n) d 


J a A ab (k)J b + J a M ab (k)U b , 


(6.9) 


where A and B are regular matrices of order N. J a means J a (—k ) and, in matrix notation, 
J a can be understood as a row matrix. For example, in our case, the effective action (5.1) 
reads 


where 


J“ = 


S(2) = 


v 2 r du 

2 ./ 2tt 


j“A„ t (u.) J b + 


(a^\ 

4 ? 

\x (0 7 


( a l 1 ) 


n“ = 


Vx {3) , 


A = 


/ 0 -n 
—// —2mo 0 
^ 0 0 0 ; 



( 6 . 10 ) 


( 6 . 11 ) 


With (6.8) the action (6.9) becomes 


0 ( 2 ) _ 1 
^ren ^ 


d d k r 


41 


w>o (2 tt 
d d k 


d 


r [A ab (k) + M ac O c i (S- 1 )i(k)} J b 


(27T. 


i d 


tcl/^R jb 
J Lr ab J 


( 6 . 12 ) 


where the range of u is chosen to be positive following the prescription in [13]. Notice that 
©“(S -1 )], is independent of J“, because neither §“ nor ©“ depends on J a . Furthermore 
©“(§ _1 )k is independent of the choice of the initial condition (6.3), because the different 
choice of initial value vectors are nothing but a linear transformation tpf —> which 

induces right multiplications in the solutions: § —> § R, O —> OR. This proves the Green 
functions are independent of choice of basis for our initial conditions. 

Notice that since A and B are also independent of J, G bi b is independent of J and man¬ 
ifestly gauge invariant, giving alternative reason for the invariance of the Green functions 
under the residual gauge symmetry. 


7 Conclusion 

We investigated the gauge invariance of physical observables in a holographic theory under 
the local diffeomorphism. We find that gauge invariance is closely related to the holo¬ 
graphic renormalisation. Apart from the zero-th mode residual gauge dependence, gauge 
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dependence is cancelled by the local counter terms defined in the boundary, However, due 
to the difference in the space-time structure between the near-horizon and near boundary 
regions, there are residual gauge structure near boundary. There is a subtle and deep 
connection between the degrees of freedom at the boundary and those at the bulk. There 
are three degrees of freedom at the boundary, out of which only two can be embedded into 
bulk fields such that they are the boundary values of the bulk fields satisfying the incoming 
boundary conditions. The residual gauge symmetry is not a proper gauge symmetry but a 
solution generator near the boundary. We proved the invariance of Green’s functions under 
such a symmetry in the context of algorithm by which all AC transports are constructed 
simultaneously. 

We can extend the RGS such that it satisfies the infalling boundary condition, which 
we call the boudanry shifting diffeomorphism. Then we can view things more concisely and 
natural. RGS is not gauge symmetry but a solution generating transformation. Therefore 
it generate formally new solution both in boundary and bulk. By extending it to BSD, the 
bulk part of the solution can be accepted as a true bulk degree of freedom more naturally 
since the latter satisfies the in-falling boundary condition. 

Note added: After this work is almost finished, the paper [14] appeared where residual 
gauge invariance was discussed using a different method. 
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